Isotope-selective photo-ionization for calcium ion trapping 
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We present studies of resonance-enhanced photo-ionization for isotope-selective loading of Ca + 
into a Paul trap. The 4s 2 1 S *-* 4s4p 1 ~P 1 transition of neutral calcium is driven by a 423 nm laser 
and the atoms are photo-ionized by a second laser at 389 nm. Isotope-selectivity is achieved by 
using crossed atomic and laser beams to reduce the Doppler width significantly below the isotope 
shifts in the 423 nm transition. The loading rate of ions into the trap is studied under a range 
of experimental parameters for the abundant isotope 40 Ca + . Using the fluorescence of the atomic 
beam at 423 nm as a measure of the Ca number density, we estimate a lower limit for the absolute 
photo-ionization cross-section. We achieve loading and laser-cooling of all the naturally occurring 
isotopes, without the need for enriched sources. Laser-heating/cooling is observed to enhance the 
isotope-selectivity. In the case of the rare species 43 Ca + and 46 Ca + , which have not previously been 
laser-cooled, the loading is not fully isotope-selective but we show that pure crystals of 43 Ca + may 
nevertheless be obtained. We find that for loading 40 Ca + the 389 nm laser may be replaced by an 
incoherent source. 

PACS numbers: 32.80.Fb, 32.80.Rv, 32.80.Pj 



I. INTRODUCTION 



Resonance-enhanced photo-ionization for loading ion 
traps was first demonstrated with magnesium and cal- 
cium ions by Kjaergaard et al. Q. The same group has 
recently used crossed-beam Doppler-free excitation of the 
272 nm 4s 21 S <-*■ 4s5p 1 P i calcium transition, followed 
by photo-ionization, to load selectively all the naturally 
occurring isotopes of this element; the loading was not 
observed directly, but by using charge-exchange to in- 
fer the presence of the other isotopes from the fluores- 
cence of laser-cooled 40 Ca + ions which replaced them 
Guide et al. showed that, using the 4s 2 1 S <-> 4s4p 1 P 1 Ca 
transition at 423 nm, followed by excitation close to the 
continuum by an ultraviolet photon, photo-ionization is 
around five orders of magnitude more efficient than con- 
ventional electron bombardment ionization [j| . 

As pointed out by these authors, photo-ionization has a 
number of advantages over electron bombardment. Only 
the desired species is loaded into the ion trap, allow- 
ing pure crystals of particular isotopes to be obtained. 
Since no electron beam is involved, there is no charging 
of insulating parts of the trap structure, which leads to 
drifting electric fields that for many experiments must 
be accurately compensated. The efficiency of the photo- 
ionization process allows much lower number densities of 
neutral atoms in the interaction region, greatly reducing 
the quantities of material sputtered onto the trap elec- 
trodes; it has been shown that clean electrode surfaces 
reduce the heating rate of trapped ions from the motional 
ground state Q , an important consideration for quantum 
logic experiments in ion traps. For experiments in which 
high-finesse optical cavities are combined with ion traps, 
the lower atomic beam density reduces the degradation of 



the mirror surfaces and can avoid the need for a separate 
"loading trap" 0,0. 

Our particular interest lies in using 43 Ca + as an ion- 
qubit in quantum logic experiments iTl. This ion has 
a number of advantages as a qubit 0, Oj ; it is also an 
attractive candidate for a trapped ion optical frequency 
standard |ToL Iit| . The most obvious difficulty in work- 
ing with 43 Ca + is its low natural abundance of 0.135%. 
In this paper we describe the development of a photo- 
ionization system capable of loading 43 Ca + into a Paul 
trap from a natural abundance source. (Isotopically- 
cnriched sources are available, although the maximum 
enrichment of ~ 80% means that it would still be ad- 
vantageous to use an isotope-selective method for reliable 
loading of this isotope.) 

We describe the crossed atomic and laser beam setup 
for excitation of the 423 nm 4s 2 1 S <-* 4s4p 1 P 1 transi- 
tion with a Doppler width narrow compared with the 
isotope shifts. Spectroscopy of the atomic beam is used 
to ensure the intensity of the 423 nm laser is below sat- 
uration (to avoid reducing the isotope-selectivity of this 
step by power broadening), and also to estimate the num- 
ber density of atoms in the beam. Photo-ionization from 
the 4s4p 1 P 1 level is achieved by a 389 nm photon; since 
we expect the photo-ionization cross-section to be only 
weakly dependent on the photon energy, we use a non- 
stabilized diode laser for this step. We study the photo- 
ionization trap loading rate as a function of the power 
and detuning of the 423 nm laser, and as a function of the 
power of the 389 nm laser. We find that the 389 nm laser 
photo-ionizes equivalently above or below lasing thresh- 
old, and that it can be replaced by a much cheaper in- 
coherent source if high photo-ionization rates are not re- 
quired. We estimate the absolute photo-ionization load- 
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ing efficiency and compare it with that for electron bom- 
bardment ionization in our apparatus. 

The sensitivity of the loading rate as a function of 
the 423 nm detuning determines the isotope-selectivity 
which can be attained. However, for very rare species 
with mass number X a practical limitation is the rate 
at which charge-exchange replaces the desired, trapped, 
ions X C& + by the abundant 40 Ca + ions issuing from the 
oven. We demonstrate loading of all the naturally occur- 
ring isotopes, including 43 Ca + and the 0.004%-abundant 
46 Ca + . These two species are laser-cooled for the first 
time in the experiments reported here. The limited 
isotope-selectivity attainable for these rare ions means 
that some form of "purification" of the ions loaded is 
necessary in the case of 43 Ca + we show that this 

is possible with no loss of the desired species. Hence it 
should be possible to load arbitrarily large crystals of this 
important ion. 

Level diagrams for 43 Ca and 43 Ca + , showing the en- 
ergy levels and transitions relevant to these experiments, 
are shown in figure ^ Isotope shifts in the six transi- 
tions shown are summarized in table [IJ together with the 
abundances of the six naturally occurring isotopes. 

II. APPARATUS 

A. Ion trap 

The apparatus consists of a linear radio- frequency (r.f., 
6.2 MHz) Paul trap in an ultra-high vacuum system 
(background pressure < 2 x 10 -11 torr), with diode laser 
systems available to drive all the transitions shown in 
figure n Details of the trap are given in [3. Typical 
axial and radial frequencies used in the present exper- 
iments are 170 kHz and 500 kHz respectively. Auxiliary 
d.c. electrodes are used to compensate stray electric fields 
to the level of about 1 V/m using a photon-r.f. correla- 
tion technique applied to the fluorescence from a single 
trapped ion [5l| • A magnetic field of ~ 3 G is applied to 
prevent optical pumping into magnetic sub-states 

B. Lasers 

The lasers are all of the external-cavity grating- 
stabilized Littrow design, except the 389 nm laser which 
has no grating stabilization. The blue and ultravio- 
let laser diodes are all GaN devices, the infrared ones 
GaAlAs. The grating-stabilized lasers may be locked to 
tunable, stabilized, optical cavities to reduce linewidths 
and medium-term frequency drift below 5 MHz. The 
laser frequencies are set within about 100 MHz of the 
ionic transitions with wavemeters. Doppler-free satu- 
rated absorption spectroscopy in a calcium hollow cath- 
ode is used to monitor the 423 nm laser frequency rel- 
ative to the 4s 21 S <-> 4s4p 1 P 1 atomic transition. The 
vacuum wavelength of the 389 nm free-running laser was 



measured using a grating spectrograph to be 388.9(1) nm 
at 20° C (with negligible dependence on drive current); 
above lasing threshold, the spectral width is expected 
to be < 100 MHz, below threshold ~3nm. For con- 
venience of alignment, the 423 nm and 389 nm laser 
beams are both injected into the same single-mode op- 
tical fibre to transport them to the trap; this ensures 
that the beams are well overlapped in the interaction 
region, which is centred on the trap. The transmis- 
sion of the fibre is ~ 30% at 423 nm and ~ 10% at 
389 nm. The maximum intensities at the centre of these 
beams are ~50mW/mm 2 at 423 nm and ~5mW/mm 2 
at 389 nm. Typical intensities of the cooling lasers at 
the trap are: 10mW/mm 2 at 397 nm; 20mW/mm 2 at 
393 nm; 6mW/mm 2 at 866 nm; 200mW/mm 2 at 854 nm; 
650mW/mm 2 at 850 nm. These intensities are based on 
the measured spot sizes for the various beams, which 
range from 40-300 /im. 



C. Detection system 

Fluorescence from the interaction region at 423 nm and 
397 nm acts as a diagnostic of neutral calcium and cal- 
cium ions respectively. The trap region is imaged by a 
compound lens onto an aperture to reject scattered light; 
further lenses re-image the trap, via a violet filter, onto a 
photomultiplier (PMT). A movable beamsplitter allows 
a portion of the light to be directed to a charge-coupled 
device (CCD) camera; the camera is used to detect the 
presence of non-fluorescing ions in the trap (revealed as 
"gaps" in an otherwise regular crystal), and to check 
whether the ions form a cold crystal or a hot cloud. The 
net detection efficiency of the PMT system, including 
the solid angle subtended by the lens and the quantum 
efficiency of the PMT, was measured to be 0.13(2)% at 
397 nm; the photon count rate observed from a single, 
cold, 40 Ca + ion under conditions of near-saturation is 
consistent with this figure. The net detection efficiency 
at 423 nm is rj = 0.15(3)%. 

The discriminated pulses from the PMT are counted 
directly by 10 MHz on-board counters in the computer 
which controls the experiment; two alternately-gated 
counters are used to eliminate read-out delay. For the 
large peak photon counting rates encountered in studies 
of the 423 nm fluorescence (up to ~ 10 6 s _1 ), we correct 
the measured PMT count rate S' for the maximum count- 
ing rate S n , ax — 10 MHz, assuming S <C S n , ax , to obtain 
the true count rate S » S'/{1 - S"/S' max ). 

D. Calcium beam 

The calcium source consists of an oven made from 
stainless steel tube (diameter 2 mm, wall thickness 
0.1mm) filled with granules of calcium metal, closed by 
crimping at each end, and with an orifice (estimated area 
~lmm 2 ) filed in the side facing towards the trap. The 




FIG. 1: Energy levels of interest in 43 Ca (for photo-ionization) and 43 Ca + (for laser-cooling), with hyperfine levels labelled by 
the total angular momentum F. The nuclear spin is I =_?■ The largest hyperfine splitting in each term is indicated in MHz; 
the levels approximately conform to an interval rule |15lll6t|. Diagrams for the even isotopes are the same but with no hyperfine 
structure, (a) 43 Ca: the lifetime of the 4s4p 1 P i level is 4.49(4) ns [r^| and the natural linewidth of the 423 nm transition 
35.4(3) MHz. Photo-ionization from the 4s4p 1 P i level requires a photon of wavelength less than 389.8 nm. (b) 43 Ca + : the 
natural widths of the ultraviolet transitions are approximately 23 MHz, the lifetimes of the metastable D levels about 1.2 s [l8|. 
Inset: Positions of the hyperfine components of the 393 nm and 397 nm transitions (solid lines, labelled with F lovcr : F uppor ) 
and the isotope shifts of the even isotopes (dotted lines, labelled with the mass number), as a function of detuning A relative 
to 40 Ca + . The isotope shifts in the infrared transitions (similar in all three transitions) are also shown, without the hyperfine 
details for 43 Ca + ; note the change of scale. 



oven is heated by passing a current (between 3 A and 
6 A) along its length, which produces a beam of calcium 
atoms directed towards the trap. Estimated oven tem- 
peratures are in the range 500-650 K, depending on the 
current used (section llll C|) . 

The oven is situated 22(2) mm from the centre of the 
trap. The atomic beam is collimated by the r.f. trap elec- 
trodes to an effective aperture of length / = 1.7(1) mm 
along the direction of the 423nm/389nm laser beams. 
However, the detection system only images a limited 
length of the interaction region (inset figure |2J, so the 
effective aperture for observation of the 423 nm fluores- 
cence is reduced to I' = 0.44(4) mm. Since the oven 
orifice is of comparable dimensions to I, V the angular 
divergence of the beam is somewhat larger than these 



distances would otherwise imply. 



III. SPECTROSCOPY OF THE ATOMIC BEAM 

We describe in this section spectroscopy of the 423 nm 
transition, using the calcium atomic beam from the oven. 
For optimal isotope-selectivity in the photo-ionization 
loading, the narrowest possible linewidth of this transi- 
tion is required. Ideally the linewidth should be small 
compared with the isotope shifts (table QJ, with the 
lower limit being set by the 35.4(3) MHz natural width. 
Other contributions to the homogeneous width include 
power broadening (saturation), transit-time broadening 
and laser linewidth. Collisional broadening is entirely 
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Mass 


Natural 


Isotope shifts (MHz) 


number 


abundance 


CcL So Pi 


Ca+ S-P 


Ca+ D-P 






423 nm 


397 nm [393 nm] 


854 nm [850, 866 nm] 


40 


96.9% 











42 


0.647% 


394 


425(6) 


-2350(4) 


43 (e.g.) 


0.135% 


612 


688(17) 


-3465(4) 


44 


2.09% 


774 


842(3) 


-4495(4) 


46 


0.004% 


1160 


1287(4) 


-6478(8) 


48 


0.187% 


1513 


1696(6) 


-8288(7) 



TABLE I: Abundances of the naturally occurring calcium isotopes |T^| . and isotope shifts for Ca and Ca + transitions relevant 
to this work. Shifts in the 423 nm transition are taken from [Tfij ]: uncertainties quoted are below 1 MHz. Shifts in the ultraviolet 
Ca + transitions were measured in |20|l : we quote the (more accurate) measurements for the 397 nm transition: the shifts at 
393 nm are all consistent with these. Shifts in the infrared transitions are from |15| ; the data given are for the 854 nm transition, 
which are again consistent with those for the other two lines. For the odd isotope, the shift of the centre of gravity of the 
hyperfine components of each transition is given; note that in the Ca + S-P lines, the hyperfine structure is significantly larger 
than the isotope shift, figure [T]^b). 



negligible at the low number densities involved here. 
The inhomogeneous broadening (Doppler width) is de- 
termined by the oven temperature, the collimation angle 
of the atomic beam and the angle between the laser beam 
and the atomic beam. 

For orthogonal beams, the velocity distribution ap- 
proximates to a Gaussian whose width compared with 
that in the oven is determined by the collimation geom- 
etry and the oven orifice |22j, §10.1]. For non-orthogonal 
beams, the velocity distribution broadens and becomes 
asymmetric, resulting in a broadening and shift when it 
is convolved with the homogeneous profile. 



A. Alignment of the crossed beams 

To ensure orthogonal alignment of the laser and atomic 
beams, we minimize the width of the 423 nm fluorescence 
spectrum as a function of the angle between the beams. 
The power of the 423 nm laser P423 is set safely below the 
saturation level for maximum sensitivity (see next sec- 
tion). A fluorescence spectrum with the alignment opti- 
mized is shown in figure^] together with a simultaneously 
acquired Doppler-free hollow cathode spectrum. From 
several such spectra we find that the offset between the 
40 Ca fluorescence peak and the hollow cathode saturated- 
absorption peak is at most 10 MHz (after accounting for a 
pressure red-shift of up to 13 MHz in the hollow cathode 
lamp due to 6 torr of neon buffer gas [1^]). This indicates 
a beam misalignment of at most ~ 12mrad. 

The fitted Voigt profile shown has Lorentzian full- 
width at half- maximum (fwhm) 39 MHz and Gaussian 
fwhm 52 MHz. The transit-time broadening is estimated 
from the measured laser spot size and the approximate 
atomic velocity (section flll C|l to be 2(1) MHz. The laser 
linewidth estimated from the error signal when the laser 
is locked to a stabilized reference cavity is 4(1) MHz. The 
Gaussian width is consistent with the collimation aper- 



ture I' = 0.44(4) mm and the estimated oven temperature 
T~ 610 K if we model the oven as an extended source of 
size ~ 1 mm. 

B. Saturation broadening 

It is important to avoid saturation broadening of the 
423 nm transition. Since it is difficult to estimate ab- 
solute intensities accurately, and the atomic response is 
averaged over the laser beam profile, we studied the sat- 
uration of the transition empirically, by varying the laser 
power. Results, obtained after alignment of the crossed 
beams, are shown in figure |21 Saturation effects be- 
gin, both in the peak fluorescence and the Lorentzian 
linewidth, above a laser power of about 10 fiW. 

The photo-ionization rate will increase as the popula- 
tion of the 4s4p 1 P 1 state, i.e., linearly until saturation 
sets in. Since the atomic beam is much larger than the 
laser beam, however, the total number of atoms in the 
upper state will continue to increase even well above sat- 
uration, though more slowly with increasing power. The 
isotope-selectivity will decrease roughly as the square of 
the Lorentzian linewidth, for fixed Gaussian width (the 
wings of the Voigt profile are determined overwhelmingly 
by the Lorentzian component, for the Gaussian width 
involved here). To maximize the photo-ionization rate 
without sacrificing isotope-selectivity, we work at powers 
in the range 5-15 fxW for photo-ionization loading. 

Specifically, the fractional population of the 
4s4p 1 P 1 Mj = upper state as a function of de- 
tuning A423 at A = 423 nm is, for linearly-polarized laser 
light of intensity /, 

-IAB 

where A is the Einstein coefficient for the transition 
(A = 2.23 x 10 s s -1 ; the weak decay route to 4s3d 1 D 2 
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FIG. 2: (a) Fluorescence spectrum from the atomic beam 
at an oven current of 5.0 A and a laser power below satura- 
tion, Pa23^ 6 /^W. The laser was locked to a tunable, stabi- 
lized, optical cavity to minimize drift during the scan; the 
frequency scale was determined using a second, confocal, cav- 
ity (600 MHz free spectral range). The best- fit Voigt pro- 
file is shown (solid line), where the amplitude, linear base- 
line, frequency offset, Lorentzian and Gaussian widths were 
all floated simultaneously. The total fwhm is 76 MHz. The 
abundances, isotope shifts and hyperfine splittings were given 
the known values [Trl llfij. Right inset: Magnification of the 
region containing the three hyperfine components of 43 Ca, 
whose expected positions are indicated. The dashed line 
shows the best fit when the 43 Ca components are omitted 
from the theoretical profile. The largest contribution to the 
signal by 43 Ca is about 22% of the total signal above back- 
ground. Left inset: Fluorescence at 40 Ca line centre from 
the interaction region, as observed on the CCD camera. The 
imaging system aperture limits the length of the 423 nm beam 
which is observed to 0.44(4) mm. (b) Simultaneous Doppler- 
free saturated absorption signal from the hollow cathode 
lamp; logarithmic transmission is shown. The Doppler-free 
peak (fwhm ~ f 70 MHz) is superimposed on the Doppler- 
broadened background (fwhm ~ 2.2 GHz, dashed curve), and 
is centred within 10 MHz of the 40 Ca fluorescence peak. Peaks 
due to the other isotopes are too small to be visible. 
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FIG. 3: Fluorescence and Lorentzian linewidth on the 423 nm 
transition, as a function of laser power P423, obtained from 
fitting 423 nm spectra with Voigt profiles (whose Gaussian 
width was fixed at 51 MHz). The oven current was 4.25 A. 
The vertical dashed lines indicate the power at which the in- 
tensity at the centre of a beam with 100 /im spot size reaches 
I sei = 2%hcA/\ 3 = 3.68mW/mm 2 . (a) Peak (corrected) 
PMT count rate above background on the 40 Ca component. 
A straight line (y oc x) is fitted to the data for P423 < 10 fiW; 
the onset of saturation is visible as the signal drops below 
this line at high powers, (b) Analysed Lorentzian fwhm for 
the same data. The linewidth starts to increase for power 
above about 10 /^W. The horizontal dashed line indicates the 
35 MHz natural linewidth. The curve shows the expected ho- 
mogeneous fwhm, obtained by averaging the atomic response 
over a circular Gaussian beam with 100 fim spot size and as- 
suming a 4 MHz laser linewidth. 



r L = 4 MHz, is plotted in figure EJb). 



C. Number density and temperature in the atomic 

beam 



may be neglected since the calculated lifetime against this 
decay is long compared to the transit time across 
the laser beam), (3 — (27iT L + A)/2 with I\ the laser 
linewidth, the intensity is measured in units of a satura- 
tion intensity J sat = 2whcA/\ 3 and A 42 3,r L are in units 
of Hz. This expression must be integrated over the spa- 
tial profile of the laser beam to obtain the expected signal 
from the atomic beam; the resulting homogeneous fwhm, 
assuming a circular Gaussian beam and a laser linewidth 



We can use the peak °Ca fluorescence at 423 nm to de- 
duce the number density of calcium atoms in the atomic 
beam, using the known line strength and the measured 
photon collection efficiency. This in turn permits an esti- 
mate of the absolute photo-ionization loading efficiency. 
By measuring the number density as a function of oven 
current we can also compare the efficiency of photo- 
ionization trap loading with the electron bombardment 
method (which have similar loading rates at very differ- 
ent oven currents). 
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The peak fluorescence as a function of oven current 
is shown in figure Ela). The laser power was set below 
saturation and under these conditions the peak corrected 
photon count rate S is related to the number density TV 
of atoms in the ground state by 



S = 



3X 3 A 
8irhc 



A23/V • V(0) • a 4a N 



where 040 = 96.9% is the 40 Ca abundance, I' = 
0.44(4) mm the length of the interaction region observed 
through the imaging aperture (measured along the laser 
beam direction), 77 = 0.15(3)% the net collection effi- 
ciency of the detection system at A = 423 nm, r — 1.43(2) 
a geometrical factor arising from the angular distribu- 
tion of the fluorescence and V(0) the value of the fitted 
Voigt profile at line centre (which is normalized such that 
J °° V(y — vq)&v = 1). The fitted values of V(0) show a 
slight decrease with increasing oven current, due to in- 
creasing inhomogeneous width as the temperature rises, 
but this is no greater than the ~ 10% fitting uncertainty, 
and for the purpose of figure EJa) we use the average 
value V(0) = 9.9(1.0) GHz -1 so that the number density 
7V[m~ 3 ] = 1.3(3) x 10 6 S'[s" 1 ] and it can be plotted on 
the same ordinate. 

For a small interaction region on axis of the atomic 
beam, the number density N in the interaction region is 
related to that in the oven Nq by 



N = 



Npa 
4nd 2 



where a is the area of the oven orifice and d the distance 
from the orifice to the interaction region [24l §4.2]. In 
our apparatus, d — 22(2) mm but a is not accurately 
known. However, since number density depends strongly 
on temperature T, we can make a reasonable estimate 
of T by assuming that a — 1.0(5) mm 2 , the value sug- 
gested by our profile analysis. We then find a temper- 
ature which satisfies p — NokT and the known vapour 
pressure curve p — p(T) simultaneously [25j . This is 
shown in figure Efb) > an d ranges between about 500 K 
and 650 K. The temperature is usefu l for estim ating the 
mean velocity along the beam v = y/ irkT/2M for atoms 
of mass M:v = 433(4) m/s for 40 Ca at T = 575(10) K. 
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FIG. 4: Calcium number density and oven temperature as 
a function of oven current, (a) Corrected photon count rate 
above background on the 40 Ca peak obtained from fitting 
423 nm spectra with Voigt profiles, where the Lorentzian 
width was fixed at 41 MHz from a spectrum at 4.25 A oven 
current. The right-hand ordinate gives the calcium number 
density in the interaction region deduced from this fluores- 
cence rate as described in the text, and the error bars refer 
to the number density. A quadratic regression fit is shown 
(solid line). At 3 A, the signal was so small that the line 
shape was poorly determined and the Gaussian width was 
also fixed for the fit. (b) Estimated oven temperature based 
on the number density in (a), with quadratic regression (solid 
line). The dashed lines indicate the temperature range im- 
plied by the 50% uncertainty assigned to the oven orifice area 
a = 1.0(5) mm 2 . 



quantitative estimate of this is difficult since the angular 
divergence depends on the sizes of the oven orifice and the 
trap capture region, neither of which is accurately known. 
Instead, we investigate the actual photo-ionization load- 
ing rate as a function of 423 nm detuning, as described 
in the following section. 



D. Conclusion 



With the 423 nm beam well aligned orthogonal to the 
atomic beam, and the laser power set below saturation, 
we observe a linewidth of 76 MHz (fwhm) , made up of 
Lorentzian and Gaussian contributions of about 39 MHz 
and 52 MHz respectively. This compares favourably with 
the isotope shifts, table |U We do not use these data to 
predict relative loading probabilities for the different iso- 
topes, however, since the collimation angle of the atomic 
beam for loading the trap is greater than that relevant to 
the 423 nm fluorescence (I > I') and the Gaussian contri- 
bution to the linewidth will be correspondingly larger. A 



IV. PHOTO-IONIZATION LOADING STUDIES 

A. Expected form of photo-ionization cross-section 

Calculations [2g of the cross-section Q for photo- 
ionization from the 4s4p state suggest that, at 389 nm, 
Q lies in the range 60-280 megabarns (lmegabarn = 
10~ 22 m 2 ). The large range results from the existence of 
an auto-ionizing resonance close to the ionization thresh- 
old, whose wavelength is uncertain to ~ 5 nm. Since the 
energy of the resonance is not well-known, we do not at- 
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tempt to choose a particular laser wavelength, beyond 
ensuring that it is below the threshold required for ion- 
ization. The width of the resonance is certainly large 
compared to the laser linewidth, so there is no need for 
the laser to be carefully frequency-stabilized. 

The electric fields present in a Paul trap can allow ion- 
ization to occur even for excitation slightly below the 
continuum limit 3]. For the typical velocity estimated 
above, atoms cross the photo-ionization laser beams in 
~ 0.2 /is and all atoms are thus able to sample the ex- 
trema of the r.f. field (period 0.16/is). The peak field 
varies across the trap, but is at most ~ 10 5 V/m. A 
rough estimate [22j of the effect of such a field on the 
ionization threshold may be made by assuming a hydro- 
genic potential V(x) — —e/Aiteox near the continuum 
and equating this to the potential (—eEx) of the elec- 
tron in the external field E; the result is that the ion- 
ization limit could be shifted at the furthest to a wave- 
length of 391.3 nm. However, no significant enhancement 
of the photo-ionization cross-section is expected by using 
wavelengths closer to this shifted ionization limit, even 
for monochromatic light, because the oscillatory electric 
field blurs out the discrete energy levels of the atom near 
the continuum and the average transition probability is 
no higher than that within the continuum |23, §18.6]. 



B. Photo-ionization loading rate measurements 

A typical photo-ionization loading curve is shown in 
figure El where the fluorescence at 397 nm from cold 
40 Ca + ions in the trap is plotted as a function of time 
after the photo-ionizing lasers are switched on. As the 
size of the ion crystal grows, the fluorescence per ion 
increases due to the effect of r.f. micromotion for ions ly- 
ing off-axis, because they spend time closer to resonance 
with the red-detuned 397 nm cooling laser. However, for 
the range of crystal sizes in these experiments, the fluo- 
rescence per ion is roughly constant and we use the final 
fluorescence level as a measure of the total number of ions 
loaded. (The fluorescence was calibrated after a typical 
load by reducing the trap strength until the ions could 
be counted on the camera.) For each load, we wait a few 
seconds after the photo-ionizing lasers are switched off, 
to allow any hot ions to be cooled and to join the crystal, 
check that the ions are still crystalline using the cam- 
era, and then take the final number of ions m divided by 
the duration r of exposure to the photo-ionizing lasers 
as our measure of the loading rate R. We find that an 
error s/m/r is broadly consistent with variations in the 
loading rate taken under nominally identical conditions 
(see for example figure EJ). 

We have studied the loading rate as a function of the 
power P423 and detuning A423 of the 423 nm laser, and 
of the power P389 of the 389 nm laser. This was done 
at fixed oven current, to keep the atomic number den- 
sity as constant as possible. For each measurement, the 
397 nm fluorescence was allowed to increase to a level rep- 
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time (s) 

FIG. 5: Typical photo-ionization loading fluorescence curve. 
The 423 nm and 389 nm photo-ionizing lasers are switched 
on at t = and off at t = 216s (dotted line). The 397 nm 
and 866 nm cooling lasers illuminate the ions continuously. 
Background counts due to scattered 397 nm light and (dur- 
ing the loading) 423 nm fluorescence and 389 nm light have 
been subtracted from the PMT signal. The step in 397 nm 
fluorescence when the first cold ion appears is arrowed. The 
large steps in the fluorescence are typical and may result from 
the details of the cooling dynamics as the number of ions in 
the trap increases. The dashed line indicates the final fluo- 
rescence level, which corresponds to 36(3) ions; the average 
loading rate of 0.17(3) ion/s is obtained by dividing this by the 
time for which the photo-ionizing lasers are switched on. The 
conditions were: oven current 4.25 A, Pi23~ 6 ^tW, A423 = 
0(5) MHz, P 38 9 = 0.93(9) fxW, A397 = -68(10) MHz. 



resenting about 40 trapped ions, whereupon the photo- 
ionizing lasers were blocked. In practice, the number of 
ions loaded varied between 22(2) and 63(5); the length 
of photo-ionizing exposure was between 2 s and 600 s. 

Figure E] shows the dependence of loading rate on 
423 nm power and detuning. Note that the loading rate 
could in principle saturate at a different power level to the 
423 nm fluorescence studied above, since photo-ionization 
is limited by the width of the 389 nm beam (the lifetime 
of the 4s4p 1 P 1 level, 4.49(4) ns, is much shorter than 
the time an atom spends in the laser beams, ^0.2^s) 
whereas the fluorescence collected is limited by the imag- 
ing aperture. Nevertheless, we observe that the loading 
rate also starts to saturate above a power of P423~ 10 fiW 
due to saturation of the upper state population. Given 
the power available in the 423 nm laser, the maximum 
loading rate could be made perhaps an order of mag- 
nitude larger by increasing the spot size at the trap to 
1mm; spot sizes much larger than this, however, would 
start to reduce the isotope-selectivity because of the in- 
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FIG. 6: (a) Photo-ionization loading rate as a function of 
423 nm laser power, with: oven current 4.25 A, A423 = 
0(10) MHz, P 3 89 = 0.93(9) mW,A 397 = -80(20) MHz. The 
loading rate starts to saturate for P423 > 10 fiW. The curve 
is drawn to guide the eye. (b) Loading rate versus 423 nm 
detuning, with: oven current 4.25 A, Pi23~ 6 /jW, P389 = 
80(8) /iW, A397 = -97(10) MHz. A Voigt profile (solid line) is 
fitted to the data; the frequency offset, amplitude and Gaus- 
sian width were floated whilst the Lorentzian width was fixed 
at 39 MHz from the fit in figure The fwhm of the fitted 
curve is 94 MHz. The frequency origin is taken from the fit- 
ted offset of a 423 nm fluorescence spectrum taken during the 
same run; the fitted profile here shows a blue-shift of 12 MHz 
which probably results from the slightly different collimation 
geometry relevant to loading. 



creasing atomic beam divergence angle. 

The loading rate as a function of A423 is of critical 
importance for isotope-selective loading. We expect the 
dependence to be modelled by a Voigt profile with sim- 
ilar Lorentzian width to that of the 423 nm fluorescence 
spectrum, but different Gaussian width because of the 
different atomic beam collimation angles for fluorescence 
detection and ion loading. The Lorentzian and Gaussian 
widths are not well-determined if they are both floated, 
due to the noise in the data, so we fix the Lorentzian 
width at 39 MHz, from the fluorescence spectrum in fig- 
ure |3 and find a fitted Gaussian width of 71 MHz. This 
Gaussian width is consistent with a beam collimation an- 
gle defined by the electrodes and an oven orifice of size 
~ 1 mm. The consequent isotope-selectivity expected is 
discussed in section IV~Al 



FIG. 7: Photo-ionization loading rate as a function of 389 nm 
laser power, with: oven current 4.25 A, Pt23~ 6 /iW, A423 = 
0(5)MHz,A 3 97 = -68(10) MHz. No significant difference in 
efficiency is observed even when the laser is well below lasing 
threshold (vertical dashed line), indicating that the photo- 
ionization cross-section is independent of the spectral width 
of the laser. A straight line (y oc x) is fitted to the data. The 
single measurement with the ultraviolet LED in place of the 
389 nm laser is also included (see section HV PI . 



In contrast to the 423 nm bound-bound transition, no 
saturation effect is expected for the 389 nm bound-free 
ionizing transition, and this is verified for the 389 nm 
powers we can attain (figure 0): the dependence of the 
loading rate on P389 is linear over two orders of mag- 
nitude. Significantly, there is no change in efficiency 
whether the 389 nm laser is running above or below las- 
ing threshold, despite the fact that the spectral width 
is expected to be about five orders of magnitude larger 
(~ 3 nm) below threshold. This implies a weak depen- 
dence on wavelength of the photo-ionization cross-section 
in this region. 



C. Absolute photo-ionization efficiency 

We can deduce a lower limit for the photo-ionization 
cross-section from the loading rate under conditions of 
near-saturation on the 423 nm transition, when we ex- 
pect that nearly half the atoms are in the 4s4p 1 P 1 upper 
level (at least across that part of the 423 nm beam which 
overlaps the 389 nm beam). From the definition of the 
photo-ionization cross-section Q, the number of atoms 
ionized per unit time is 

R = Ql ■ \a i0 N ■ P 389 /e 

where I = 1.7(1) mm is the length of the interaction re- 
gion defined by the trap electrodes, (^a4oiV) the num- 
ber density of 40 Ca atoms in the 4s4p 1 P 1 level, P389 the 
389 nm laser power, and e the 389 nm photon energy. 

To minimize uncertainty in the number density due 
to variations in oven conditions, TV was determined 
shortly before measuring the loading rate and found to 
be 1.7(4) x 10 11 m -3 . At maximum 423 nm power, the 
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estimated intensity at the centre of the laser beam was 
14(4)J sat which, for |A 423 | < 10 MHz, gives n 4a4p > 
0.48(1) within the half-intensity radius of the beam 
(equation justifying the assumption of saturation. 
The loading rate was measured to be R — 4.9(9) ion/s 
with P389 = 1.1(1) //W, giving a cross-section Q = 
170(60) megabarns. This is a lower limit since it assumes 
that every ion produced is trapped and that the capture 
region of the trap extends across the full width of the 
atomic beam as far as the r.f. electrodes, but it lies in 
the anticipated range (section llV A|) . 

Alternatively, the loading rate can be expressed in 
terms of the probability q per atom of being ionized (and 
trapped) as it crosses the interaction region: 

R = q ■ ^a i0 N ■ lw 2 /t 

where w 2 is the cross-sectional area of the 389 nm beam 
(which limits the extent of the interaction region), and 
t s» w/v is the time an atom of velocity v spends in the 
interaction region. For the conditions above and the laser 
beam size w ~ 0.1 mm, we have t 0.2 /is and find q 
8 x 10 -7 . The maximum q we observed (with maximum 
powers in both the photo-ionizing lasers) was q « 4 x 
10 -5 , and by eliminating the optical fibre and using the 
maximum available 389 nm power ~2mW, we would in 
principle be able to obtain a maximum q rs 10~ 3 , i.e., 
one in 1000 atoms crossing the interaction region would 
be loaded. 



D. Loading using an ultraviolet LED 

The observation that the 389 nm laser is able to photo- 
ionize atoms just as efficiently below lasing threshold led 
us to experiment with an alternative light source for the 
second step of the photo-ionization, an ultraviolet light- 
emitting diode (LED). The device tested was a Nichia 
NSHU590 diode, with nominal output power 500 fiW, 
peak wavelength 375 nm and spectral width 12 nm. The 
chief advantages of this device are that it is about three 
orders of magnitude cheaper than a wavelength-selected 
ultraviolet laser diode, and does not require temperature 
and current stabilization. Measurements with a grat- 
ing spectrograph indicated that the peak wavelength was 
in fact between 375 nm and 390 nm, depending which 
part of the emission region was focussed on the spectro- 
graph slit, and the spectral width was ^25nm; however 
at least half of the spectral intensity lies in the range 
A < 389.8 nm useful for photo-ionization. The main 
drawback of the LED is its extended emission region: 
it is difficult to focus a large fraction of the power into 
a region smaller than a few millimetres, unlike a laser 
beam. 

When we first tested the LED in place of the 389 nm 
laser we were surprised to find that it was capable of load- 
ing ions into the trap while the 423 nm laser was blocked 
and the oven was switched off! A small number of 40 Ca + 
ions were loaded, together with several "dark" ions; at 



the same time we noticed that the stray d.c. electric field 
in the trap changed by ~ 10%, drifting gradually back 
to its former value over a period of an hour or so. Our 
interpretation is that the large (~ 5 mm) patch of light 
from the LED was irradiating the r.f. electrodes, ablat- 
ing and photo-ionizing calcium and other ions from the 
surface where they had previously been deposited by the 
atomic beam. Photo-electric emission could be respon- 
sible for charging insulating patches on the electrodes, 
causing the stray field to change temporarily before the 
charge leaks away (the photo-electric cut-off wavelength 
for calcium is 433 nm). Surface ablation is a well-known 
method of loading ion traps [2!j, though it is perhaps 
surprising at the low light intensities involved here. 

To avoid this problem, we imaged the LED onto a 
200 y«m pinhole, which was re-imaged at approximately 
1:1 magnification into the trap region, and overlapped 
with the 423 nm beam. This eliminated loading of 
non- 40 Ca + ions and effects on the stray field, but only 
0.32(3) /zWof ultraviolet power was then available at 
the trap. The loading rate under these conditions, 
0.06(1) ion/s, is plotted in figure and is similar to that 
expected with comparable 389 nm laser power. This rate 
is not sufficient for loading the rarer isotopes, where it 
is necessary to beat losses due to charge- exchange (sec- 
tion EU, but is perfectly adequate for loading small 
crystals of 40 Ca + . We note that high-power (~ 150 mW) 
ultraviolet LEDs are also available, though their cost ap- 
proaches that of a laser diode. 



V. ISOTOPE-SELECTIVE LOADING 

A. Expected selectivity 

From the data in figure EIb) we expect good isotope- 
selectivity, because the net fwhm, 94 MHz, is small com- 
pared with the isotope shifts in the 423 nm transition. 
This discrimination has to be set against the very low 
abundance of some isotopes. In figure [8] we plot the cal- 
culated probability of loading the different isotopes as 
a function of detuning A423, based on the fitted Voigt 
profile in figure Efb). We regard the loading as "isotope- 
selective" for those isotopes where there exists a detun- 
ing at which one is more likely to load the desired iso- 
tope than any other; it can be seen from figure |S] that 
this is the case for all isotopes except 43 Ca + and 46 Ca + . 
The most demanding case is 46 Ca + , where the maximum 
achievable loading probability, relative to the total load- 
ing probability, is only 6.0%. The relative probability of 
loading 43 Ca + is also plotted in the figure: it reaches a 
maximum of 21% at a detuning of 0.623 GHz. This com- 
pares with a theoretical maximum of 33% which would 
be possible with a Lorentzian width equal to the natural 
linewidth and negligible Gaussian width; thus not a great 
deal would be gained by efforts to collimate further the 
atomic beam. 

The detunings of the Doppler-cooling lasers also affect 
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FIG. 8: Probability of loading different calcium isotopes as 
a function of detuning of the 423 nm laser, normalized to 1 
for 40 Ca + at zero detuning. The curves are generated using 
Voigt profiles with the same width parameters as the fit in 
figure |(Jb) (except that the Gaussian width is adjusted to 
take account of the different isotopic masses). The dashed 
curves for the even isotopes, and the solid curve for 43 Ca + (the 
three hyperfine components have been summed), are labelled 
by the mass numbers. The thick solid curve is the sum of 
the contributions from all the isotopes. The detached solid 
curve gives the probability of loading 43 Ca + relative to the 
total loading probability at the same detuning, plotted in the 
region where this quantity is > 0.05. 



the selectivity, due to the isotope shifts in the ionic transi- 
tions (tableHJ. For example, to cool 42 Ca + efficiently, the 
397 nm laser needs to be detuned ~ 350 MHz to the blue 
of the 40 Ca + transition, which will tend to heat 40 Ca + 
ions and expel them from the trap, thus enhancing the 
selectivity. 

B. Charge-exchange loss 

It is possible for an ion x Ca. + which has been loaded 
into the trap to undergo charge- exchange with a neutral 
atom in the atomic beam; since the overwhelming prob- 
ability is that the neutral atom is 40 Ca the process is 
usually: 

x Ca+ + 40 Ca^ x Ca + 40 Ca+ 

and the original ion in the trap is replaced by 40 Ca + . 
The charge-exchange rate is proportional to the number 
density in the atomic beam. This process is a limitation 
for loading more than one ion of a rare isotope: if the loss 
rate through charge-exchange exceeds the loading rate, 



then rare ions already trapped are replaced by 40 Ca + 
faster than subsequent ones can be loaded. If the charge- 
exchange rate is sufficiently fast, then rare ions will not 
even be cooled before they are replaced, and the trap will 
fill up with 40 Ca + before any fluorescence signal from the 
rare ions is seen. 

In practice we find that, providing one is not attempt- 
ing to load large crystals, charge- exchange is only a limi- 
tation for the rarest isotope 46 Ca + . For the typical load- 
ing parameters used for loading 43 Ca + , for example, the 
loading rate is ^lion/s while the lifetime against loss 
due to charge- exchange is ~ 1 min. 46 Ca + is about thirty 
times less abundant, so we expect the loss rate to ap- 
proach the loading rate under the same conditions. 

We note that the charge- exchange process was used 
by Mortensen et al. ,2j to demonstrate isotope-selective 
loading, but that they were not attempting to retain and 
to laser-cool the different isotopes; instead they relied on 
charge-exchange to replace every isotope by 40 Ca + and 
used the resulting 40 Ca + fluorescence as a measure of the 
original load size. 



C. Even isotopes 

Loading one of the even isotopes consists of choosing 
the 423 nm detuning for maximum selectivity and set- 
ting the appropriate detunings of the 397 nm and 866 nm 
lasers used for Doppler-cooling. After loading, we can 
check for the presence of other isotopes by looking for flu- 
orescence using different detunings of the cooling lasers 
(the lasers are blocked while the detunings are changed, 
to avoid possible heating) . Phase-sensitive detection with 
a chopped 866 nm beam can be used to detect very small 
fluorescence signals above scattered 397 nm background 
light. 

As expected from figure|S|we can load small, pure crys- 
tals of 44 Ca + and 48 Ca + (and of course 40 Ca + ) without 
difficulty. Predicted peak loading probabilities of 44 Ca + 
and 48 Ca + are 94% and 85% respectively. Practical load- 
ing rates can be achieved in spite of the small abundances 
by using maximum 389 nm power (~ 100 /zW) and by in- 
creasing the oven current. We concentrate here on the 
more challenging cases of 42 Ca + (which lies closest to 
40 Ca + ) and 46 Ca + (which is extremely rare). 

The maximum expected probability of loading 42 Ca + 
is 58%, at a detuning A 423 = 0.394 GHz, with 40 Ca+ by 
far the dominant impurity. However, when we loaded 
using this detuning, we obtained a crystal consisting of 
~ 19 42 Ca + ions and only ~ 5 40 Ca + ions, i.e., a loading 
efficiency of about 80%. We attribute this to the mech- 
anism mentioned above, that the 397 nm laser was blue- 
detuned from the 40 Ca + transition by about 350 MHz, 
thus tending to heat this species and expel it from the 
trap; the 40 Ca + ions are only cooled sympathetically by 
the 42 Ca + j2jj. This mixed-species crystal also illustrates 
the typical charge- exchange rates: after ~ 20 min, with 
the photo-ionizing lasers blocked but the oven left on at 
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4.25 A, the number of 40 Ca + ions had increased to ~ 18 
and there was only a very low signal visible at the 42 Ca + 
detunings; after a further 30 min, there was no detectable 
signal from 42 Ca + . This indicates that the lifetime (per 
ion) against loss by charge-exchange is ~ 17 min, for a 
number density in the atomic beam of 9(2) x 10 10 m~ 3 . 
The limited isotope-selectivity achievable could be over- 
come by selective heating methods [13J, [l4| to obtain a 
pure crystal of 42 Ca + . 

The 0.004%-abundant 46 Ca+ is the most difficult iso- 
tope to load, with an expected peak loading probabil- 
ity of 6.0% (neglecting enhancement by heating as ob- 
served for 42 Ca + ). To achieve the maximum photo- 
ionization rate without excessive saturation-broadening 
of the 423 nm transition, we set ^423^ 180 //W; we used 
maximum 389 nm power (~100/iW) and increased the 
oven current to 5.0 A. Using the measured loading rates 
for 40 Ca and the number density data we expect a 46 Ca + 
loading rate of ~ 0.1 ion/s, with a sacrifice of a factor ~ 3 
in selectivity because of the high 423 nm power. Under 
these conditions, a fluorescence signal was detected af- 
ter ^17s of exposure to the photo-ionizing lasers. The 
oven was switched off immediately to prevent charge- 
exchange. The ions could not be made to crystallize, 
however, indicating that too many impurity ions were 
present to be cooled sympathetically by the 46 Ca + . We 
therefore applied a "tickle" voltage to one of the trap 
end-caps, close to the axial resonance frequency, which is 
less likely to expel the directly-cooled 46 Ca + than other 
isotopes. Crystallization was then observed, figure^ and 
camera images implied that the crystal consisted of a sin- 
gle 46 Ca + ion and several impurity ions (most probably 
40 Ca + ). The observed isotope shift in the 866 nm transi- 
tion was —6.4(2) GHz; although this particular shift has 
not been measured before, it is in agreement with that 
expected (table HJ and thus confirms the identity of the 
ion. 



D. Odd isotope 

The laser requirements for Doppler-cooling 43 Ca + are 
more demanding than for the even isotopes because of the 
hyperfine structure, figure^b). In particular, a separate 
rcpumper is necessary on one of the ultraviolet transi- 
tions because of the large (3.2 GHz) hyperfine splitting 
between the S 3 1/2 and S*, 2 levels (where the superscript 
denotes F), 

We choose the 393 nm S* /2 «-> Pg, 2 transition for cool- 
ing, both because it is an intrinsically strong F <-> F + 1 
transition and (thanks to the inverted hyperfine struc- 
ture) it is blue-detuned from the S 1/2 *-* P 3 / 2 transitions 
in all the even isotopes (inset figure This second fact 
means that any other isotopes that are loaded will tend to 
be heated. Furthermore, this transition is nearly-cycling 
(apart from the ~ 6% decay route to the D levels and 
off-resonant driving of the 160MHz-detuned S* /2 <-> Pg, 2 
component), so that the second ultraviolet laser needed 
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FIG. 9: Fluoresence signal from a sing le 46 Ca+ ion, trapped 
in the company of several "dark" ions (probably 40 Ca + ), as 
a function of 397 nm detuning relative to 40 Ca + . Three scans 
were averaged to produce the plot. The high-frequency "edge" 
in the fluorescence curve, characteristic of scanning past the 
resonant frequency, has been assigned the known isotope shift 
of 1287 MHz. The first dip in the fluorescence (arrowed) is the 
point at which the ions crystallize and the sing le 46 Ca+ ion 
becomes visible on the camera (inset). The second dip, at 
~ 1200 MHz, is probably due to a two-photon dark resonance 
between the 866 nm and 397 nm transitions |3l| . The crys- 
tallization dip, and the movement of the sing le 46 Ca+ ion 
between discrete positions in the trap after heating and re- 
crystallization, indicate the presence of other ions. 



for repumping out of the S 3 1/2 level can still repump ef- 
ficiently at low intensities: this is an advantage since it 
is necessarily red-detuned relative to the even isotopes 
and can thus cool them. We choose to use the 397 nm 
SJ /2 «-» P* /2 transition for this repumping. 

Repumping from the D 3/2 and D 5/2 manifolds is ac- 
complished with the 850 nm and 854 nm lasers respec- 
tively: the high intensities available in these beams mean 
that they can repump adequately in spite of the hyperfine 
structure which spans ~ 400 MHz in each transition. 

To load we set the appropriate detunings of these 
four lasers relative to 40 Ca + , and set A423 « 600 MHz. 
The 423 nm power is typically ~10//W, to avoid 
saturation-broadening, and the 389 nm power maximized 
(~ 100 /«W) for fastest loading. Under these conditions 
and with an oven current of 5.0 A, the loading rate is 
~1 ion/s and the lifetime against loss through charge- 
exchange ~ 1 min. This gives time to detect a single ion 
of Ca + and switch the oven off before it is replaced by 
one of 40 Ca + . With this procedure, we can reliably load 
single ions of 43 Ca + and only rarely find that we have 
loaded 40 Ca + as well; this indicates that we achieve sig- 
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nificant enhancement over the estimated maximum load- 
ing probability of only 21%, presumably because 40 Ca + 
ions are heated by the 393 nm laser and expelled from 
the trap. The maximum fluorescence rate detected from 
a single 43 Ca + ion was 23000 s -1 , compared with typ- 
ically 32000 s -1 for 40 Ca + (a reduction attributable to 
less efficient repumping from the D states). 

To confirm the identity of the first single ion we loaded 
by this method, we scanned the frequency of the 397 nm 
repumper: at low intensity, the two hyperfine compo- 
nents Sf /2 <-» P* /2 and S 3 1/2 <-» P a 1/2 are well-resolved, fig- 
ure an d have the expected separation. As a further 
check, the axial resonance frequency was measured and 
compared with that of 40 Ca + , giving a mass number of 
42.8(2). 

If the photo-ionizing lasers and oven are left on for 
longer we find that we load crystals containing about 50% 
43 Ca + and 50% 40 Ca + . This may occur either because of 
charge-exchange, or because once 43 Ca + is trapped and 
cooled, it can cool 40 Ca + sympathetically in spite of the 
heating effect of the 393 nm laser. We discovered that it 
was possible to purify such mixed crystals by temporar- 
ily blocking the 397 nm ultraviolet repumper; when this 
laser is blocked a small fluorescence signal is still visi- 
ble from 43 Ca + (because of off-resonant repumping on 
S?/ 2 P3/2 4 by the 3GHz-detuned 393 nm laser) but 
there is no longer any laser capable of cooling 40 Ca + di- 
rectly. We assist the 40 Ca + heating by switching on the 
866 nm repumper, tuned to be resonant with the 40 Ca + 
transition. One might expect the Coulomb coupling be- 
tween the different species to make this purification pro- 
cess rather inefficient, but we found on the contrary that 
it always seems successful at removing the 40 Ca + ions 
without loss of any 43 Ca + . For example, a mixed crystal 
consisting of three 43 Ca + ions and eight 40 Ca + ions was 
purified completely by blocking the 397 nm repumper for 
~ 8 s. The technique can also be applied during loading: 
with the oven and photo-ionizing lasers on, we watch the 
43 Ca + ions being loaded on the camera and, as soon as 
there is evidence of "dark ions" in the trap, block the 
397 nm repumper for a few seconds to purify the crystal, 
then un-block it to continue loading. In this way, pure 
43 Ca + crystals of arbitrary size can be loaded; the inset 
of figure HUI shows a 9-ion linear crystal. 



VI. PHOTO-IONIZATION AND ELECTRON 
BOMBARDMENT COMPARED 

The limitations of electron bombardment ionization 
are discussed in A significant problem in our ap- 
paratus was the effect of firing the electron gun on the 
stray electric fields in the trap. After loading, the static 
field necessary to compensate the micromotion of the ions 
was found to decay exponentially with a time constant 
~ 6 hours. For the typical load shown in figure lllf al 
the vertical stray field would not have been stable to 
the 1 V/m resolution of our compensation method until 
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FIG. 10: Fluorescence signal from a single 43 Ca + ion 
as the frequency of the 397 nm ultraviolet repumping 
laser is scanned. The 397 nm intensity was reduced to 
~0.2mW/mm 2 for this scan, to reduce the widths of the 
peaks. A double Lorentzian is fitted to the data: the sepa- 
ration between the two components is 582(6) MHz (where a 
1% error is allowed in deducing the frequency scale from a 
confocal cavity of 300 MHz free spectral range) , in excellent 
agreement with the known hyperfine splitting in the P1/2 level 
of 581 MHz. The S?/ 2 «-> P1/2 component has been assigned 
the known shift of -1381MHz relative to 40 Ca+. Inset: A 
pure crystal of nine 43 Ca + ions, loaded using periodic block- 
ing of the 397 nm repumper to remove even-isotope impurity 
ions, as described in the text. The separation between the 
closest ions is 13(1) /an. 



~ 28 hours after the trap was loaded. We ascribe this 
drift to charge on insulating parts of the trap structure 
which decays with time, and have not observed this sys- 
tematic decay after photo-ionization loading. On longer 
timescales, between loads, we observed that the vertical 
stray field changed by up to ~ 200 V/m; this is compared 
with recent photo-ionization compensation data in fig- 
ure E^b). The cause of these long-term changes in the 
stray field could have been changes in the distribution of 
material deposited on the trap electrodes when the oven 
or electron gun was fired: calcium deposited on the stain- 
less steel r.f. electrodes could give rise to static fields of 
this magnitude, due to the difference in work functions 
between these two metals. 

This highlights a further issue: based on the work 
of the NIST Ion Storage group, we expect that mate- 
rial deposited onto the trap electrodes from the oven 
and/or electron gun will contribute significantly to mo- 
tional heating of trapped ions in the sub-Doppler cool- 
ing regime relevant to quantum logic experiments |"3.l82|. 
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The estimates of calcium number density and temper- 
ature in the atomic beam given in section 1111 CI above, 
imply a flux of ~ 4x 10 9 atom/s/mm 2 on a surface normal 
to the atomic beam, at the 6 A oven current which was 
necessary for loading the trap using electron bombard- 
ment. For the same loading rate, this can be reduced 
by over four orders of magnitude using photo-ionization. 
Furthermore, the same photo-ionization loading rate can 
be maintained even when using a well-collimated atomic 
beam designed to prevent deposition of atoms on the trap 
structure. 

With maximum power in both photo-ionizing lasers, 
we measured a loading rate of 1-6(2) ion/s at an oven cur- 
rent of 3.25 A. This compares with an electron bombard- 
ment loading rate of ~ 1 ion/s at an oven current of 6.0 A. 
Extrapolating the number density data of figure 0{a) to 
6.0 A oven current indicates that photo-ionization is a 
factor of ~ 10 times more efficient than our electron 
bombardment ionization setup. A further two orders 
of magnitude enhancement of the loading efficiency in 
this apparatus would be available by (i) increasing the 
spot sizes of the photo-ionizing lasers to reduce satura- 
tion effects on the 423 nm transition and (ii) eliminating 
the optical fibre in the 389 nm beam. We note, how- 
ever, that the per atom efficiency of photo-ionization trap 
loading can be approached by techniques of surface ion- 
ization |33| . 

The photo-ionization method overcomes the limita- 
tions of electron bombardment ionization: the species- 
and isotope-selectivity ensure that unwanted ions are not 
trapped and furthermore allow studies of rare isotopes 
without enriched sources; there is negligible effect due to 
the loading on the stray fields (at the level of ~ 1 V/m 
in this trap); it works well irrespective of radial trap 
strength Q, so that compensation of stray fields for dif- 
ferent trap conditions is unnecessary; the high efficiency 
of photo-ionization means that the flux of the atomic 
beam can be attenuated by many orders of magnitude 
which reduces immensely problems associated with de- 
position of material on the trap structure. These advan- 
tages must be weighed against the complexities of two 
additional lasers, compared with an electron gun inside 
the vacuum system. 
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FIG. 11: (a) Drift of the stray field versus time after load- 
ing with electron bombardment ionization. The magnitude 
of the vertical field required to compensate horizontal ion mi- 
cromotion is plotted. The micromotion compensation here 
was performed by minimizing the width of the Si/ 2 «-» P1/2 
spectrum; we assign an error of ±3 V/m. An exponential 
decay function is fitted to the data (solid line) ; this has time- 
constant 6.3(7) hours and amplitude 88(4) V/m. No such drift 
is visible after loading by photo-ionization, at the 1 V/m level, 
(b) Variation of the steady-state stray field versus time, over 
the course of nearly a year for electron bombardment load- 
ing (squares, left) and over the course of 45 days for photo- 
ionization loading (circles, right). Each point is taken from a 
different load. For the electron bombardment data the error 
bar represents the estimated additional drift possible after the 
final compensation measurement. For the photo-ionization 
data, the estimated error in compensating the stray field is 
1 V/m; the standard deviation of the data set is 2.8 V/m, indi- 
cating some residual systematic variation (which may be due 
to the remaining atomic beam flux, or mechanical variations 
in the trap geometry over time). 



VII. CONCLUSION 



We have presented quantitative studies of photo- 
ionization ion-trap loading for calcium ions, using the 
fluorescence from the neutral atoms to estimate the ab- 
solute efficiency of this process. We find that, in our 
trap, this loading method is around four orders of mag- 
nitude more efficient than conventional electron impact 
ionization, comparably to Guide et al. Q, and does not 
give rise to stray electric fields. We have shown that the 
setup of Guide et al. may be further simplified by elim- 
inating frequency-stabilization for the 389 nm laser and, 
if high loading rates are not required, by replacing this 



laser with a light-emitting diode. We have optimized the 
photo-ionization procedure for the purposes of isotope- 
selection, demonstrating trapping of all the naturally oc- 
curring isotopes from a natural abundance source. The 
rare species 43 Ca + and 46 Ca + have been laser-cooled for 
the first time. 

The odd isotope 43 Ca + is of particular interest for ex- 
periments in quantum information processing and optical 
frequency standards, and we have demonstrated the abil- 
ity to load pure crystals of this ion. We note, however, 
that an enriched source would still be advantageous for 
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work with this isotope, since it would reduce any remain- 
ing deposition of material from the atomic beam on the 
trap electrodes to the minimum possible. An interesting 
feature of photo-ionization compared with electron bom- 
bardment ionization is that the laser beams allow it to be 
spatially-selective: this could be useful for loading spe- 
cific traps in proposed quantum computing architectures 
involving microtrap arrays |34j . The spatial selectivity 
could be extended to three dimensions by crossing the 
423 nm and 389 nm beams. The isotope-selectivity would 
also be invaluable for loading "refrigerator" ions |3{| in 
schemes where continuous cooling is provided by one iso- 
tope, while quantum logic operations are carried out us- 
ing another isotope: following the estimate in [2^] , the de- 
coherence rate, due to off-resonant photon scattering, of 
a 43 Ca + qubit being cooled sympathetically by a 42 Ca + 
refrigerator ion is ~ 0.1 /s for a trap with a heating rate 
of 10 3 motional quanta per second p|. 

Finally, in comparison with the work of Mortensen et 
al., it is clear that the 272 nm 4s 2 1 S <-* 4s5p 1 P i transi- 
tion is much more favourable for isotope-selective ion- 
ization of calcium than the 423 nm transition used in 
this work, since it possesses a narrower natural width, 
8.5(5) MHz ^t|, an d isotope shifts about a factor of two 
larger However, this scheme requires a significantly 
more complex laser source than the blue diode laser used 
here, and is somewhat less efficient since it involves ex- 



citation deeper into the continuum. An unfortunate cir- 
cumstance at the time of writing is the unavailability of 
further diodes at 423 nm from Nichia Corporation but 
we are optimistic that, with blue laser diodes under de- 
velopment by several other manufacturers and research 
laboratories, the situation will improve. 
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